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REVIEW  OF  RECENT  DEVELOPMENTS  IN  EVAPORATION  DUCTING  ASSESSMENT 

by 
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NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152-5000 
USA 


SUMMARY 

Continued  interest  in  accurately  assessing  and  exploiting  evaporation  ducting 
effects  has  resulted  in  some  recent  work  summarized  in  this  paper. 

Available  evaporation  duct  height  and  path  loss  models  are  compared  (xef— irf C 
This  includes  both  sensitivities  of  duct  height  values  to  meteorological  inputs  and 
agreements  between  different  path  loss  models.  The  models  compared  are  those 
.developed  by  Jeske  Hitney  ^(xef-3')^''  and  Rotheram  Reasonable  agree- 

ment of  duct  height  calculations  and  path  loss  values  derived  from  the  different 
models  was  obtained  when  applied  to  statistically  averaged  meteorological  inputs. 
Commonly  used  point  measurements  were  found  inadequate  for  reliable  evaporation  duct 
assessment.  This  led  to  the  suggestion  of  a new  shipboard  deployable  sensor  using  a 
single  temperature  sensor  for  air-sea  temperature  .measurements.  . 


-Both  long  time  statistical  duct  height  averages  and  individual  data  sets  often 
contain  unusually  high  duct  heights  not  supported  by  propagation  measurements.  V'Paulus 
(ref  5)  analyzed  the  meteorological  conditions  and  found  a bias  in  commonly  used 
measurement  techniques  toward  stable  conditions.  He  proposed  a correction  algorithm 
and  demonstrated  its  utility  convincingly  on  a set  of  radiometeorological  measure- 
ments. \ 


The  optimum  frequency  for  utilization  of  the  evaporation  duct  for  long  range 
surface  detection  appears  to  be  in  the  10-2C  GHz  band.  This  is  illustrated  by  cal- 
culating surface  detection  ranges  for  hypothetical  radars  equivalent  in  performance 
parameters  but  operating  at  different  frequencies  (ref  6} . 


PREFACE 

The  vertical  decrease  of  relative  humidity  from  100%  at  the  air-sea  boundary 
results  in  a rapid  vertical  decrease  of  the  refractive  index.  Electromagnetic  waves 
traveling  close  to  the  water  surface  may  be  significantly  influenced  by  such  a refrac- 
tive index  distribution  and  in  extreme  cases  may  be  very  efficiently  ducted.  This 
oceanic  evaporation  duct,  which  is  found  to  some  degree  nearly  all  the  tine  over  all 
oceanic  areas,  is  of  particular  importance  for  over-the-horizon,  surface-to-surface 
radar  coverage.  It  is  an  important  factor  in  the  low  flying  ship-to-ship  missile 
problem.  To  illustrate  some  fundamental  properties  of  the  effect  of  the  oceanic 
evaporation  duct  on  radio  propagation,  two  vertical  distributions  of  the  modified 
refractivity  profile  are  shown  in  figure  1.  (Modified  refractivity  M = (n-1)  - 10  4- 
h/a  - 10°;  n = refractive  index.,  h = height  above  the  surface,  a = earth's  radius.) 
The  dashed  line  represents  no  evaporation  ducting  while  the  solid  line  is  repre- 
sentative of  good  ducting  conditions.  Figure  2 shows  the  vertical  distribution  of 
path  loss  for  the  two  refractivity  distributions  assuming  a 19-nmi,  9.6  GHz  over-water 
propagation  path  with  one  terminal  at  16  feet  above  the  water  (path  loss  is  the  ratio 
in  dB  of  transmitted  to  received  power  assuming  isotropic  antennas) . The  dashed  curve 
in  figure  2 indicates  a decrease  of  path  loss  (or  an  increase  of  signal)  with  height. 
The  solid  curve  corresponding  to  the  solid  K-curve  in  figure  1 shows  a minimum  in  path 
loss  (or  a maximum  in  received  signal)  at  a height  of  about  15  feet.  At  this  height 
the  signal  enhancement  over  no  ducting  conditions  is  62  dB.  Signal  enhancements  of 
such  a magnitude  emphasize  the  importance  of  the  oceanic  evaporation  duct.  (The  solid 
curve  in  figure  2 must  not  be  extrapolated  beyond  the  height  plotted.  It  will  not 
cross  the  "no  ducting"  curve  but,  depending  on  horizon  range,  increase  with  altitude.) 
The  presentation  of  figure  2 may  also  be  used  to  determine  antenna  heights  for  optimum 
utilization  of  ducting  conditions.  For  example,  an  antenna  at  a height  of  64  feet 
would  receive  11  dB  less  signal  than  an  antenna  at  15  feet.  Nonetheless,  the  signal 
enhancement  from  evaporation  ducting  would  be  30  dB  for  the  high  antenna.  A fact 
frequently  overlooked  is  that  signal  enhancements  from  evaporation  ducting  are  usually 
encountered  for  all  possible  antenna  heights  on  board  ship. 

A convenient  parameter  for  the  description  of  ducting  conditions  is  "duct  height1 
defined  as  the  height  at  which  the  M-curve  attains  its  minimum  value.  Thus,  the  solid 
curve  in  figure  1 has  a duct  height  of  47  feet.  In  general,  greater  duct  heights  will 
result  in  higher  signals.  However,  this  duct  height  parameter  should  not  be  misin- 
terpreted as  a height  limit  above  which  signal  enhancements  cannot  occur-  As  shown  in 
figure  2,  30  dB  signal  enhancements  are  experienced  well  above  the  duct  height. 
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REFRACTIVITV  PROFILES 


The  vertical  distribution  of  the  refractive  index  must  be  known  for  propagation 
calculations-  A direct  measurement  of  the  refractive  index  profile  under  operational 
conditions  is  impractical.  For  instance,  the  profile  in  figure  1 changes  20  M-units 
within  the  first  4 feet  above  the  water  and  only  an  additional  4 M-units  between  4 and 
47  feet.  Even  for  moderate  winds  and  ocean  waves  it  would  be  very  difficult  to  make 
measurements  that  close  to  the  water  surface.  In  addition,  profiles  based  on  single 
measurements  usually  show  strong  variations  which  often  make  the  true  profile  shape 
unrecognizable.  These  variations  are  smoothed  out  only  if  mean  values  are  derived  by 
averaging  over  several  minutes.  If  the  functional  relationship  iis  known  which  governs 
the  vertical  distribution  of  the  refractive  index,  then  only  measurements  at  the  sea 
surface  and  at  a convenient  reference  height  are  needed.  Such  measurements  are  called 
bulk  measurements.  Advances  in  boundary  layer  research  have  made  it  possible  to 
derive  profile  descriptions  which  are  adequate  to  describe  radio  propagation  condi- 
tions. Those  profiles  are  dependent  on  turbulent  transfer  processes  and  assume 
different  shapes  depending  on  the  thermal  structure  in  the  boundary  layer.  Under 
thermally  neutral  conditions  (adiabatic  lapse  rate)  a logarithmic  profile  describes 
the  vertical  distribution  of  the  meteorological  parameters  of  interest.  For  non- 
neutral temperature  conditions  (diabatic  lapse  rate)  the  vertical  gradient  of  the 
meteorological  parameters  is  modified  by  a function  which  depends  on  stability. 
According  to  JesV..  (ref  2),  in  the  stable  region  a stability  function  first  proposed 
by  Monin  d':d  Obukiiov  (ref  7)  should  be  used  which  results  in  a logarithmic-linear 
profile.  Under  intense  thermal  stratifications  the  turbulence  theories  used  are  no 
longer  valid.  Therefore,  for  stabilities  exceeding  a bulk  Richardson's  number  of  one, 
np  general  functional  relationship  for  the  vertical  distribution  of  the  meteorological 
parameters  is  presently  available.  Bulk  Richardson's  number  can  be  determined  from 
the  temperature  difference  AT  (in  K)  between  the  surface  and  a reference  height  z (in 
m)  and  the  wind  speed  u (in  knots,  measured  at  z.)  according  to  1 


Rib  " 1*3zi 


AT 


(1) 


The  logarithmic-linear  profile  also  provides  the  satisfactory  profile  descrip- 
tions in  the  unstable  region  (R..<0) . However,  a different  stability  function 
provides  more  accurate  profiles  for  superadiabatic  lapse  rates  (ref  2) . Based  on  the 
two  stability  functions  for  stable  and  unstable  conditions  duct  height  DH  (in  m)  may 
be  calculated  according  to  the  following  formulas. 


Stable  region  (0<Rib<l) : 


DH  = AN/  (b^B-AN  a/X,'  ) 


(2) 


where 


B = ln(Zl/z0)  + cz-j/L' 

AN  is  the  refractivity  difference  between  surface  and  reference  height,  z1 , b.  is  a 
constant  with  the  value  of  -0.125  n,  the  hydrodynamic  roughness  of  the  sea  z1  = \.5  x 
10  m,  a = 4.5,  and  L'  is  the  so-called  Monin-Obukhov  length.  The  Monfn-Obukhov 
length  can  be  determined  from  an  empirical  relationship  between  bulk  Richardson's 
number  and  a so  called  profile  coefficient  (ref  2). 

Unstable  Region  (R^fa<0) : 

DH  = [(bjB/AJJ)4— 4a(b1B/A  H)3/I,')-1/4 


where 


B = ln(Z;l/z0)  - $ (3) 

The  function  A may  be  found  in  ref  8.  Further  refinements  of  boundary  layer  models 
have  been  attempted  based  on  the  so-called  Businger-Dyer  stability  function  (ref  9- 


DUCT  HEIGHT  NOMOGRAMS 


The  measurements  necessary  for  dijct  height 
(1)  - (3)  are  surface  water  temperature  TS  and,  at 
axr  temperature  TA,  relative  humidity  RH,  and  wind 
rated  from 
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Where  P is  atmospheric -pressure  in  millibars,  T is  temperature  in- Kelvins,  and  e is 
water  vapor  pressure  in  millibars.  A simple  nomogram  (figure  3)  was  developed  which 
provides  an  easy  means  for  deriving  evaporation  duct  height  based  on  simple  bulk 
meteorological  measurements.  Starting  at  graph  1 in  figure  3,  a horizontal  line  is 
drawn  for  the  measured  RH  and  TA  values  until  it  intersects  the  appropriate  TS  curve 
in  graph  2.  From  this  intersection  a vertical  line  is  drawn  into  graph  4 where  the 
intersection  with  a horizontal  line  starting  at  the  measured  TA-TS  and  WS  values  of 
graph  3 gives  the  duct  height  DH.  Since  evaporation  duct  height  by  itself  does  not 
give  quantitative  systems  performance  data,  additional  nomograms  have  been  developed 
in  which  the  graph  for  duct  height  is  replaced  by,  for  instance,  a radar  detection 
range  graph.  In  figure  4 the  range  of  detecting  a 30000  m target  with  a 90%  prob- 
ability using  an  SPS-10  radar  is  hhown  (graph  4).  Besides  eliminating  cumbersome  and 
lengthy  calculations,  figures  3 and  4 illustrate  also  the  sensitivity  of  evaporation 
ducting  assessment  to  the  measurement  accuracy  of  the  input  parameters.  For  example, 
small  changes  to  air-sea  temperature  differences  values  may  significantly  change  duct 
heights  and  detection  ranges.  Both  nomograms  shown  in  figures  3 and  4 are  based  on 
models  and  their  verification  described  by  Hitney  (ref  3)  and  Richter  and  Hitney  (ref 
11). 

MODEL  SENSITIVITIES 

Patterson  (ref  1)  recently  compared  the  relative  performance  of  evaporation  duct 
and  path  loss  models.  He  found  that  the  models  performed  reasonably  well  when  con- 
sidering statistically  averaged  meteorological  inputs  but,  not  surprisingly,  were 
extremely  sensitive  to  point-observed  meteorological  input  parameters.  He  concluded 
that  the  measurement  techniques  employed  by  naval  and  transiting  commercial  vessels 
are  not  of  sufficient  quality  to  calculate  individual  evaporation  duct  heights. 
Figures  5 and  6 illustrate  the  spread  between  measured  and  calculated  9.624  GHz  path 
loss  data  based  on  the  models  used  by  Rotheram  (ref  4)  and  by  the  Integrated 
Refractive  Effects  Prediction  System  (IREPS)  (ref  2,  11,  12,  13) . (The  measured  data 
are  from  ref  11,  November  period,  path  length  19  nmi,  transmitter  antenna  height  5 m) . 
Since  small  uncertainties  in  air-sea  temperature  differences  (see  figure  3)  may  cause 
significant  changes  in  duct  height,  an  evaporation  ducting  sensor  development  program 
is  being  conducted  at  the  Naval  Ocean  Systems  Center,  which  uses  a single  sensor  for 
TA  and  TS  measurement.  This  eliminates  the  need  of  accurate  cross  calibration  for 
different  sensors.  A very  fast  response  thermistor  probe,  connected  by  wire  to  a 
shipboard  processing  unit,  is  thrown  from  a convenient  shipboard  location  into  water 
undisturbed  by  the  moving  ship.  During  the  fall  the  probe  measures  TA  and,  sub- 
sequently, TS  when  entering  the  water.  The  probe  may  be  retrieved  or  made  expendable. 
Initial  tests  showed  that  air-sea  temperature  differences  of  0.1  Celsius  are  ob- 
tainable. This,  by  far  exceeds,  commonly  achievable  accuracies  and  should  improve 
individual  evaporation  duct  height  measurements  significantly. 
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CORRECTION  ALGORITHM 

Examining  evaporation  ducting  statistics  like  the  one  shown  in  figure  7 (for 
Marsden  Square  88,  Hawaii),  one  often  finds  high  percentages  of  duct  height  occur- 
rences for  large  duct  heights  (all  duct  heights  in  excess  of  40  m are  combined  in  the 
column  40-42  m) . Such  large  evaporation  duct  heights  are  difficult  to  explain  from  a 
micrometeorological  interpretation,  are  not  supported  by  radio  propagation  observa- 
tions and  are  not  found  in  meteorological  measurements  for  which  special  care  has  been 
taken  like  meteorological  buoy  observations  (ref  5) . The  most  likely  error  source  is 
surface  heating  in  TA  measurements  which  leads  to  erroneous  high  thermal  stabilities. 
Paulus  (ref  5),  based  on  a number  of  careful  arguments,  introduced  a modified  duct 
height  calculation  algorithm  in  which  stable,  high  humidity  conditions  are  accepted 
and  stable,  low  humidity  conditions  are  assumed  to  be  in  error  and  defaulted  to  un- 
stable. This  algorithm  applied  to  the  data  of  figure  7 produces  a much  more 
reasonable  distribution  and  is  shown  in  figure  8.  Even  more  convincing  is  an  applica- 
tion of  Paulus'  (ref  5)  technique  to  18  GHz  propagation  measurements  over  an  81  km 
over-water  path  (transmitter  and  receiver  heights  were  19  m and  11  m,  respectively) . 
The  solid  curve  in  figure  9 is  the  calculated  path  loss  for  this  geometry  and  the  dots 
are  measured  values.  When  Paulus'  technique  (ref  5)  was  applied  to  the  measured 
values  of  figure  9,  a significantly  better  agreement  was  obtained  as  shown  in  figure 
10.  From  the  evidence  so  far  it  appears  that  Paulus'  modified  duct  height  calcula- 
tions are  capable  of  dramatically  improving  evaporation  ducting  statistics. 

SURFACE  DETECTION  RANGES 

Based  on  measurements  described  in  ref  11,  it  appears  that  signal  enhancements 
from  evaporation  ducting  are  largest  in  the  10-20  GHz  frequency  range.  For  higher 
frequencies  molecular  absorption,  scattering  from  the  rough  ocean  surface,  in- 
homogeneities in  the  evaporative  duct  and  scattering  from  aerosols  counteract  the 
otherwise  increased  effectiveness  of  the  wave  guiding  mechanism.  Anderson  (ref  6) 
compared  four  hypothetical  radars  with  equivalent  performance  characteristics  but 
operating  at  frequencies  of  3,  6,  10,  18  GHz.  He  used  a global  refractivity  climatol- 
ogy (ref  14)  and  propagation  models  discussed  in  ref  15  to  calculate  the  probability 
of  detecting  surface  targets  at  ranges  well  beyond  the  normal  radar  horizon.  The 
results  are  displayed  in  figure  11  and  show  both  a latitude  dependence,  i.e.,  the 
evaporation  ducting  effect  is  stronger  in  lower  latitudes,  and  the  superior  perfor- 
mance for  frequencies  in  the  10-20  GHz  frequency  band. 
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Waveguide  Calculations 
Observation  Period:  1 to  24  June  82 


Figure  9.  Comparison  of 
experimental 
measurements 
to  waveguide 
predictions. 


5 ie  15  20 

Evaporation  Duct  Height  (ml 


Haveauido  Calculations 
Observation  Period:  1 to  24  June  82 
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Evaporation  Duct  Height  <m) 
Modified  Duct  Height  for  RSTD  > -1 


Figure  10.  Comparison  of 
experimental 
measurements 
to  waveguide 
predictions  in 
which  all  data 
points  based  on 
an  air-sea 
temperature 
difference  greater 
than  -1°C  were 
subjected  to  a 
modified  duct 
height  calculation. 


Probability  of  Surface  Target  Detection 


Figure  11.  Geographic 

influence  of  the 
evaporation  duct 
on  shipboard  radar 
surface-target 
detection  capabil- 
ities. The  target 
is  well  beyond  the 
normal  horizon  and 
the  increased  de- 
tection range  is 
solely  due  to  the 
evaporation  duct. 
Predictions  are 
for  four  radar 
frequencies;  3,6, 
10,  and  18  GHz. 
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CONCLUSIONS 

Evaporation  ducting  is  an  important  propagation  enhancement  mechanism  over  the 
oceans  of  continued  scientific  interest.  Recent  efforts  have  shown  that  presently 
used  models  appear  adequate  for  statistical  assessments.  An  important  correction 
algorithm  was  developed  to  remove  erroneously  high  duct  heights.  A new  air-sea  tem- 
perature probe  is  under  development  which  is  expected  to  improve  in-situ  evaporation 
ducting  assessment.  Future  work  should  address  propagation  modeling  under  horizon- 
tally inhomogeneous  conditions,  investigation  of  sensing  techniques  for  the  three- 
dimensional  evaporation  ducting  structure  and  development  of  meteorological 
forecasting  techniques. 
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DISCUSSION 


E.Vilar 

(1)  You  presented  results  obtained  in  the  early  70s  at  several  frequencies,  X-band  in  particular,  between  two  Greek 
islands.  Could  you  tell  us  the  transhorizon  path  length? 

(2)  We  are  currently  concerned  in  Europe  with  transhorizon  interference  over-sea  paths  (part  of  Eurocop  COST 
210).  Could  you  please  comment  about  the  desirability  and  minimum  sample  of  meteorological  data  to  correlate 
our  observations  with  the  existence  of  an  evaporation  duct? 

Author’s  Reply 

(1)  Path  length  was  19  nmi  with  terminal  heights  of  5 m (see  figure  16  of  reference  15). 

(2)  The  parameters  describing  the  evaporation  duct  must  be  carefully  measured  as  described  in  the  preceding  paper. 
In  addition,  upper  air  soundings  should  be  taken  to  differentiate  between  ground-based  ducting  effects  caused  by 
elevated  refractive  layers  and  evaporation  ducting  effects. 

L.B.FeIsen 

Because  the  physical  environment  presents  a multiparameter  problem,  where  the  parameters  themselves  are  not  well 
defined,  no  progress  can  be  made  on  analytical  and  numerical  modelling  unless  some  sort  of  error  bars  can  be  put  on 
the  parameters.  Those  tolerances  then  can  guide  the  modeller  as  to  whether  attempts  at  improving  presently  available 
predictive  algorithms  are  likely  to  bear  fruit.  Do  you  feel  that  more  can  or  must  be  done  to  sharpen  the  question,  or  have 
we  reached  an  impasse? 

Author’s  Reply 

More  can  and  must  be  done  in  providing  experimental  data  to  validate  existing  models  and  guide  the  development  of 
new  ones.  The  problem  in  obtaining  good  experimental  data  is  cost  and  time.  Measurement  of  all  relevant  physical 
parameters  is  expensive  and  encountering  a sufficient  large  range  of  the  uncontrolled  variables  may  require  Jong 
measurement  periods  or  many  repetitions.  Once  acquired  and  properly  documented,  good  experimental  data  can  and 
have  been  used  for  many  decades  in  various  model  validation  purposes. 

H.Vissinga 

Can  you  elaborate  somewhat  on  the  physics  of  the  evaporation  duct  in  the  presence  of  a rough  sea  surface.  In  view  of 
the  fact  that  duct  height  and  wave  height  are  not  that  different  in  magnitude  can  an  evaporation  duct  exist  over  a rough 
sea? 

Author’s  Reply 

As  one  can  see  from  figure  3,  duct  height  increases  with  increasing  surface  winds  under  unstable  conditions.  Since 
surface  winds  are  directly  related  to  sea  surface  roughness,  one  would  not  expect  the  simultaneous  occurrence  of  rough 
seas  and  small  duct  heights.  Extremely  high  winds  will  create  a well  mixed  atmosphere  with  no  evaporation  duct.  In  our 
measurements,  sea  surface  roughness  affected  path  loss  values  only  for  frequencies  well  above  X-band. 


